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Abstract. - We present mathematically rigorous expressions for the residence-time and first- 
passage-time distributions of a periodically forced Brownian particle in a bistable potential. For 
a broad range of forcing frequencies and amplitudes, the distributions are close to periodically 
modulated exponential ones. Remarkably, the periodic modulations are governed by universal 
functions, depending on a single parameter related to the forcing period. The behaviour of the 
distributions and their moments is analysed, in particular in the low- and high-frequency limits. 
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The amplification by noise of a weak periodic signal acting on a multistable system is known 
as stochastic resonance (SR). A simple example of a system showing SR is an overdamped 
Brownian particle in a symmetric double-well potential, subjected to deterministic periodic 
forcing as well as white noise. Despite of the amplitude of the forcing being too small to 
enable the particle to switch from one potential well to the other, such transitions can be 
made possible by the additive noise. For sufhcicntly large noise intensity, depending on the 
forcing period, the transitions between potential wells can become close to periodic. This 
mechanism was originally proposed by Benzi et al. and Nicolis and Nicolis [1-3] in order to 
offer an explanation for the close-to-periodic occurence of the major Ice Ages. Since then, 
it has been observed in a large variety of physical and biological systems (for reviews see, 
e.g., [4-7]). 
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Although much progress has been made in the quantitative description of the phenomenon 
of SR, many of its aspects are not yet fully understood. Mathematically rigorous results have 
so far been limited to the regimes of exponentially slow forcing [8,9], or moderately slow 
forcing of close-to-threshold amplitude [10,11]. 

One of the measures introduced in order to quantify SR is the residence-time distribution, 
that is, the distribution of the random time spans the Brownian particle spends in each 
potential well between transitions. SR is characterized by the fact that residence times are 
more; likdy to be close to odd multiples of half the forcing period than not. The residence- 
time distribution was first studied by Eckmann and Thomas for a two- level system [12]. For 
continuous systems, it has been estimated, in the case of adiabatic forcing, by averaging the 
escape rate for the frozen potential over the distribution of jump phases [13,14]. 

For larger forcing frequencies, however, the adiabatic approximation can no longer be 
used. An alternative approach is to consider time as an additional dynamic variable, which 
yields a two-dimensional problem. In the absence of noise, the system has two stable periodic 
orbits, one oscillating around each potential well, and one unstable periodic orbit, which 
oscillates around the saddle and separates the basins of attraction of the two stable orbits. 
The residence-time distribution is closely related to the distribution of first passages of the 
stochastic process through the unstable orbit. This problem was first investigated by Graham 
and Tel [15,16] and Day [17-19], and later by Maier and Stein [20] and others {e.g., [21,22]). 

At first glance, however, this two-dimensional approach seems to produce a paradoxi- 
cal result. Indeed, it is known from the classical Wentzell-Freidlin theory [23-25] that the 
distribution of first-passage locations through a periodic orbit looks uniform on the level of 
exponential asymptotics [17]. This is due to the fact that translations along the periodic orbit 
do not contribute to the cost in terms of action functional. How can this fact be conciled with 
the quasistatic picture, which yields residence times concentrated near odd multiples of half 
the forcing period? Obviously, the answer has to lie in the subexponential behaviour of the 
distribution of transitions. 

In this Letter, we extend previous results of [19-21] to a mathematically rigorous expres- 
sion for the first-passage-time distribution up to multiplicative errors in the subexponential 
prefactor, valid for a broad range of forcing periods [26,27], from which we then deduce the 
residence-time distribution. A particularly interesting aspect of the result is that both dis- 
tributions are governed by universal periodic functions, depending only on the period of the 
unstable periodic orbit times its Lyapunov exponent. All the model-dependent properties of 
the distribution can be eliminated by a deterministic time change. 

Assumptions. - We consider one-dimensional stochastic differential equations of the form 

dxt = -^V{xt,t)dt + adWt, (1) 

where Wt is a standard Wiener process, describing white noise, and the small parameter 
a measures the noise intensity (the diS'usion constant being D = cr^/2). The double- well 
potential V{x, t) depends periodically on time, with period T. The simplest example is 

V{x, t) = ^x'^ -^x^ - A sin(wi)a;, (2) 

where the forcing has angular frequency u) = 27r/T and amplitude \A\ < ■\/4/27. 

Our results apply to a general class of T-periodic double- well potentials. We assume that 
for each fixed t, V{x,t) has two minima at Xf 2{t) and a saddle at X^{t), such that 



Xl{t) < ci < < C2 < X^{t) Vt 



(3) 
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for two constants ci,C2 (in the particular case of the potential one can take C2 — — ci = 
l/\/3). Using Poincare maps, it is then straightforward to show that in the absence of noise, 
the system has exactly three periodic orbits, one of them unstable and staying between ci 
and C2, which we denote by a;P'"'(t). We introduce the notations 

a{t) = ^^V{x^^'^{tlt) (4) 
for the curvature of the potential at a;P°''(t), and 

X^^£a{t)dt (5) 

for the Lyapunov exponent of the unstable orbit. We assume that A is of order 1, but T can 
become comparable to Kramers' time. 

Finally, we need a non-degeneracy assumption for the system, which assures that the 
action functional is minimized on a discrete set of paths, and excludes symmetries other than 
time-periodicity [26]. In particular, it should not be possible to transform the equation into an 
autonomous one by a time-periodic change of variables. In the special case of the potential 
this condition is met when A Q. In addition, we will assume that \A\ is of order 1, while 

« \A\. 

Results. - Assume the system starts at time to in a given initial point in the left-hand 
potential well. The residence-time distribution is closely related to the distribution of the 
first-passage time, that is, the (random) first time r at which xt crosses the unstable periodic 
orbit xP'^'it). 

Our main result states that the probability distribution of r is governed by the following 
function, in a sense made precise in Theorem 1 below. Let 

Pit, to) = ^QMm -\h.a\)-pLe-mt)^oito)]/XT^i.) A.a„.(t,to) (6) 

N ArK(cr) 

where we use the following notations: 

• Tk(cr) is the analogue of Kramers' time in the autonomous case; it has the form 

TK(a) = ^e^/'^^ (7) 

(T 

where V is the constant value of the quasipotential on 2;P™(t). V can be computed by a 
variational method, as the minimum of the action functional over all paths connecting 
the bottom of the left-hand potential well to a::P"(f) (see [25]). In the limit of small 
forcing amplitude, V reduces to twice the potential barrier height. The prefactor has 
order cr~^ rather than 1 [21,28], due to the fact that most paths reach x'^°'^{t) through 
a bottleneck of width a (the width would be larger if J^J were not of order 1 [29]). 

• Qxriy) is the announced universal periodic function, of period AT; it has the explicit 
expression 

oo 

QAT(y) = 2Ar A{y-k\T) with A(z) = - e'^^ exp{-- e-^-}, (8) 

k— — oo 

and thus consists of a superposition of identical asymmetric peaks, shifted by a distance 
AT. The average of Qxriy) over one period is equal to 1. 
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• 6{t) contains the model-dependent part of the distribution; it is an increasing function 
of t, satisfying e{t + T) = e{t) + XT, and is given by 

9{t) ^ const +[ a(s)ds- iln^^, (9) 
Jo 2 v{0) 

where v{t) is the unique periodic solution of the differential equation v{t) ~ 2a{t)v{t) + l. 
It is related to the variance of Eq. linearized around a;P°' (i), and has the expression 

v{t) = / exp<^ / 2a{u)du\ds. (10) 



• /trans(ii^o) accounts for the initial transient behaviour of the system; it is an increasing 
function satisfying 

<^M-^Y^^^^^]) forA(i-io)<2|lna| 

(11) 

1-Oi ^ j for A(t-<o) ^ 2|lnCT| 

where L is a constant, describing the rate at which the distribution in the left-hand 
well approaches metastable equilibrium. The transient term thus behaves roughly like 
exp{-ie-[''(*)-^(*«)l /cr2(l-e-2[9(*)-e(*o)l)}. However, /trans(i,io) can be different when 
starting with an initial distribution that is not concentrated in a single point. 

• is the normalization, which we compute below. 
The precise formulation of our result is the following: 

Theorem 1. For any initial time to, any A ^ ^/a, and all times t^ to, 

rt+A 

P{r e [i,t + A]} p{s,to)ds[l + r{a)], (12) 

where r(a) — 0{^/a). 

If it were not for the limitation on A, which is due to technical reasons, this result would 
show that the probability density of r is given by p{t, to)[l + r{a)]. We expect the remainder 
to be of order a rather than ^/a■. This result has been derived in [26] in the simplified setting 
of a piecewise quadratic potential, with explicit values for V, C and L, r{a) = a, and no 
restriction on A. A full proof for the general case will be given in [27]. 

The main idea behind the proof is that sample paths reaching a;P"'(t), say, during a time 
interval \t,t -I- A] C [nT, (n -I- 1)T], are concentrated in a neighbourhood of n deterministic 
paths, the most probable exit paths, or MPEPs. Each of these paths contributes to the 
probability (|12(l . The fcth term of the sum (O is the contribution of a MPEP remaining inside 
the left-hand well for n — k periods, and then idling along a;P°' (t) during the remaining k 
periods (extending the sum from k G {0, . . . , rt — 1} to Z only results in an error of order a). 
The special form of the sum |(SJ), involving double-exponentials, has been previously noted 
in [18,20] and [21]. 
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Fig. 1 - Plots of the first-passage-time distribution p{t, 0) (full curve) for various parameter values. 
The broken curve is proportional to the average density l|14^ , but scaled to match the peak height in 
order to guide the eye. The a;-axis comprises 8 periods on each plot; the vertical scale is not respected 
between plots. Parameter values are V = 0.5, A = 1 and (a) a = 0.4 (i.e., D — = 0.08), T — 2, 
(b) a = 0.4, T = 20, (c) a = 0.5 (i.e., D = 0.125), T = 2 and (d) a = 0.5, T = 5. 



Discussion. - Let us now analyse the expression ® in more detail. 

Taking 0{t)/X as new time variable eliminates the factor 0'{t)/X in the density. Thus 
9{t)/X can be considered as a natural parametrization of time, in which one has to measure 
the first-passage-time distribution in order to reveal its universal character. We may thus 
henceforth assume that 6{t) — Xt. 

The universal periodic function Q\t depends only on the single parameter AT. For large 
AT, it consists of well-separated asymmetric peaks, while for decreasing AT these peaks overlap 
more and more and Qxriv) becomes flatter. In fact, one can easily compute the Fourier series 
of Q\T, which reads 

QMy) - E 2-''/^^r(i + ^) e'^-^'^y^'^ . (i3) 

q&I, ^ ^ 

Since the Euler Gamma function F decreases exponentially fast as a function of the imaginary 
part of its argument, Qxriv) is close, for small AT, to a sinusoid of mean value 1 and amplitude 
exponentially small in l/2Ar. 

The remarkable fact that | lncr| enters in the argument of Qxt has been discovered, to our 
best knowledge, by Day, who termed it cycling [18,19]. It means that as a decreases, the 
peaks of the first-passage-time distribution are translated along the time-axis, proportionally 
to I ln(T|. See also [20] for an interpretation of this phenomenon in terms of MPEPs. 

The remaining, non-periodic time dependence of ^ corresponds to an averaged density. 
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and behaves roughly like 



This function grows from to almost 1 in a time of order 2| Inul/A, and then slowly decays 
on the scale of the Kramers time TK(cr). It is maximal for A(< — <o) — V /a^ . 

The first-passage-time distribution is thus in effect controlled by two parameters: the quan- 
tity AT, measuring the instability of the saddle, which determines the shape of the distribution 
within a period; and the Kramers time, which governs the decay of the average density (|14|) . 

For small T/Tx{<t), the first-passage-time distribution consists of many peaks whose height 
decreases only slowly on the time scale rK(cr). Fig.Q] shows first-passage-time distributions 
for relatively large noise intensities, in order to make the decay more apparent. Increasing the 
period for constant noise intensity has two effects (Fig.n(a) and (b)): the peaks become 
narrower relatively to the period, while their height decreases faster. When T increases 
beyond rK(cr), the distribution becomes dominated by a single peak, and one enters the 
synchronization regime, with the particle switching wells twice per period. Increasing the 
noise intensity for constant period (Fig.n(a) and (c)) also produces a faster decay of the 
peaks, while at the same time the peak's location is shifted due to the cycling phenomenon. 
Fig.n](b) and (d) correspond to the same value of T/Tk{(j), but in (d) a larger noise intensity 
is responsible for broader peaks. 

Moments of the first-passage-time distribution can easily be computed up to a correction 
stemming from r(cr) (the correction due to /trans (^) is of smaller order). Using the Fourier 
series one finds (for d{t) = Xt) 



E{t"} = ^n\T^{ar 



^^''^ (l-27rigTK(a)/r)"+i^(,^ ' AT 



[l + 0(r(a))]. (15) 



In particular, taking n = yields the normalization N. Note that limcr^o o"^ logE{r} = V, in 
accordance with the classical Wentzell-Freidlin theory [23-25] . 

Two other limits are of particular interest. For AT ^ 1, the decay properties of r(l + ix) 
imply E{r"} = n!TK((T)"[l -I- 0{e~'^^ Z'^^'^) + 0{r{a))], which is close to the moments of an 
exponential distribution with expectation Tk((t). This is natural since the periodic modulation 
becomes fiat in this limit. However, it is also true that for T <C Tk{<j), one has E{t"} — 
nl TK(cr)" [1 -I- 0(T/TK(cr)) -I- 0(r(o-))], independently of the value of AT. This is due to the 
fact that the period of modulation is short with respect to the scale of exponential decay. The 
moments of the first-passage-timc distribution can thus differ significantly from those of an 
exponential distribution only when T is not too small compared to both A~^ and Tk{<t). 

A third limit in which the first-passage-time distribution should approach an exponential 
one is the limit of vanishing forcing amplitude A. However, the expression © does not hold 
in cases where A is not large compared to cr^, because it makes use of the saddle-point method 
in the vicinity of MPEPs. An asymptotic expression for A <^ has been proposed in [13,14]. 

Finally, the residence-time distribution of the system can be deduced from the knowledge of 
the first-passage-time distribution. Assume that the Brownian particle makes a transition from 
the right-hand to the left-hand potential well, crossing the unstable orbit at time s, and, having 
visited the left-hand well, crosses the unstable orbit again at time s + t. The residence-time 
distribution (i) is obtained [12] by integrating p(s -I- s)V'(s) over one period, where V'(s) is 
the asymptotic distribution of arrival phases (i.e., times modulo T). Assuming that, as for the 
potential the two wells move half a period out of phase, ■!/'(s) — Q\t{^{s—T /2)— \ hia\)/T. 
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The Fourier series ((T^ yields, for 9{t) — Xt and up to a multiplicative error 1 + 0{r{a)), the 
residence-time distribution 



Prt(t) 



1 1 

NTK{a) 



l + 2^(-l)^ 
9^1 



r 1 



Tnq 
AT 



■ f2TTqt\ 



(16) 



where f transit) has the same behaviour as /trans 0)- The periodic part of the distribution is 
minimal in multiples of T, and maximal in odd multiples of r/2, while its nonperiodic part 
behaves like the non-periodic part of the first-passage-time density. 

Conclusion. - The most important aspect of our rigorous expression for the residence- 
time distribution is the fact that it is governed essentially by two dimensionless parameters, 
AT and T/Tk{(t)^ which can be modified independently. The ratio T/Tk{<j) between period 
and Kramers time appears in most quantitative measures of SR, which indicate an optimal 
amplification when T is close to 2Tk((t). In this regime, the probability of transitions be- 
tween potential wells becomes significant during each period. The parameter AT, by contrast, 
controls the concentration of residence times within each period. Large values of AT yield a 
sharply peaked residence-time distribution, regardless of the peak's relative height. 
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